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ABSTRACT 

We address the question of the occurrence of star formation bursts in the solar neigh- 
bourhood by using the metallicity distribution of the G dwarfs. We present a method 
to recover the star formation history using simultaneously the metallicity distribu- 
tion and the age-metallicity relation. The method associates the number of stars in a 
given metallicity interval with the corresponding time interval predicted by the age- 
metallicity relation. We take into account corrections relative to observational errors, 
cosmic scatter and scale height effects. The method is tested by simulations of the 
chemical evolution of a galaxy with irregular star formation rate, and our results show 
that it is possible to recover the star formation history fairly well. The application of 
the method to the solar neighbourhood shows evidences for at least two strong events 
of star formation: a burst 8 Gyr ago, and a lull 2-3 Gyr ago. These features confirm 
some previous studies on the star formation history in the local disk. 
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1 INTRODUCTION 

The star formation rate (hereafter SFR) is one of the main 
parameters for galactic evolution studies. Previous attempts 
to recover the SFR history in the solar vicinity yield a 
present relative birthrate 6(Tg) close to unity, where To is 
the present time, and 6(Tg) is the ratio between the present 
SFR, B(Tq), and the average past SFR, (B). This fact has 
led many authors to conclude that the SFR history in the so- 
lar neighbourhood was almost constant, or slightly decreas- 
ing over the lifetime of the Galaxy (see Scalo 1986, for a 
review) . 

However, 6(Tg) ~ 1 indicates only that B(Tg) = (B), 
and not that B(t) = (B) at any time t, an expected result 
on the assumption that the SFR has been constant. There is 
no specific information on the past behaviour of the SFR in 
the function 6(Tg). Particularly, there is no indication that 
the SFR must be a smooth monotonic function of time, a 
hypothesis generally adopted for its computational simplic- 
ity or by the theoretical faith in the power of self-regulation 
in complex systems (Noh & Scalo 1990). 

In fact, several evidences strongly suggest that the local 
SFR has been quite irregular, with well-marked periods of 
enhancement and quiescence extending for about 1-3 Gyr. 
These evidences come from several independent studies rel- 
ative to (i) stellar age distributions (Twarog 1980, revisited 
by Noh & Scalo 1990; Meusinger 1991a); (ii) Ca II emis- 
sion in late-type dwarfs (Barry 1988; Soderblom, Duncan 
& Johnson 1991); (iii) distribution of Li abundances (Scalo 
1995, private communication); (iv) present-day mass func- 



tion (Scalo 1987), (v) stellar kinematics (Gomez et al. 1990; 
Marsakov, Shevelev & Suchkov 1990), and (vi) white dwarf 
luminosity function (Noh & Scalo 1990; see also Diaz-Pinto 
et al. 1994). 

The major results emerging from these studies were 
summarized by Majewski (1993). He recognizes the occur- 
rence of at least three star formation bursts, named as burst 
A, B and C. Burst A is a present star formation burst, while 
bursts B and C occurred approximately 5-6 Gyr and 8 Gyr 
ago, respectively. Between bursts A and B there is a long era 
of quiescence in the star formation activity, which is proba- 
bly reflected in the Vaughan-Preston gap in the distribution 
of chromospheric emission among late- type dwarfs (Vaughan 
& Preston 1980; Barry 1988; Henry et al. 1996). 

The occurrence of bursts is likely to leave signatures in 
the metallicity distribution of G dwarfs, in the sense that 
there will be more stars with the mean metallicity corre- 
sponding to the time of the burst. In the present work, 
we make an effort to recover the SFR history from these 
possible signatures in the metallicity distribution of nearby 
long-lived stars, as recently derived by us (Rocha-Pinto & 
Maciel 1996). In section 2, we present the method to re- 
cover the SFR history. Section 3 describes the simulations 
we have performed in order to test our method. In section 4, 
we present the results from the application of the method 
to the solar neighbourhood, and our main conclusions. 
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2 A METHOD TO RECOVER THE SFR 

We want to relate the number of stars born in the time 
interval (t, t + dt) with the total number of stars observed 
at a time T, which have metallicities between (z,z + dz). 
The first function corresponds to the absolute SFR, B(t); 
the second corresponds to the metallicity distribution of the 
region studied at time T, denoted as N(z\T). 

Let Z(t) be the mean gas metallicity, which defines an 
age-met allicity relation (AMR) for the system under study. 
The number of stars born in the interval (t,t + dt) with 
metallicities in (z, z + dz) is 



d 2 n(z, t) = B(t)P(z, t) dt dz, 



(1) 



where P(z, t) is a probability function for a star born at t to 
have a metallicity z. According to Tinsley (1975) and Basu 
& Rana (1992), 



P(z,t) = 



1 



exp 



Z(t)f 



2a 2 



(2) 



where a is the cosmic abundance scatter in the interstellar 
medium at birth. Integrating eq. (1) over time we have, 



N{z\T) dz = 



B(t) exp 



[z~Z{t)f 
2a 2 



dtdz. (3) 



It can be seen from eq. (3) that Bit) can be obtained if 
we have the metallicity distribution N(z\T) and the age- 
metallicity relation Zit). 

Assuming that the AMR is a monotonous function of 
time, and that there is no cosmic scatter in the interstel- 
lar medium, each star born at t would have a well defined 
metallicity Zit), so that we can write for a time interval At 



B(t)At = N(Z)AZ, 



(4a) 



where AZ is the corresponding metallicity variation. Eq. 
(4a) becomes for At — > 0, 



Bit)dt = NiZ)dZ. 



(4b) 



This procedure is analogous to assign a chemical age to each 
star, defined as the age for which the mean gas metallicity is 
Z. The recovery of the SFR history in our method is based 
on a distribution of chemical ages. 

However, there is a real cosmic scatter in the interstel- 
lar medium and, as a consequence, the metallicity z of a 
star is generally different from the mean gas metallicity Z 
at birth. In order to improve the distribution of chemical 
ages, we need to transform the observed metallicity distribu- 
tion N(z\T) in a distribution of mean metallicities at birth, 
N(Z\T). Hereafter, for the sake of simplicity, we will omit 
the epoch T relative to the metallicity distribution, which 
will be noted simply as N(z). Substitution of eq. (4b) into 
eq.(3) yields 



N(z) = 



rZ(T) 

I 

JZ(0) 



N(Z) exp 



2a 2 



dZ, 



(5) 



where we have also omitted the time dependence of the 
metallicity. Eq. (5) shows that N(z) is the convolution of 
the desired function N(Z) with a cosmic scatter gaussian 
function. 

Eq. (5) is much simpler than eq. (3). The method to 
solve it is the same used by Rocha-Pinto & Maciel (1996) 



to correct the observed metallicity distribution by obser- 
vational errors and cosmic scatter, and was first used by 
Pagel & Patchett (1975). The observed metallicity distri- 
bution must be approximated by a gaussian function for 
deconvolution purposes. The mean metallicity distribution 
is obtained by 

N(Z)= Niz)+8N a iz)\ z=z , (6) 
where 

5N {z) = Gliz)-Gliz), (7) 

and Go is the gaussian function adjusted to iV(«), and G? 
is its deconvolution, assuming a — 0.20 dex. 

That formalism assumes that the metallicity distribu- 
tion is complete and representative of the chemical evolution 
of the system. To eliminate some possible bias from scale 
height effects, we must consider the Sommer-Larsen (1991) 
/-factors. In this case, the mean metallicity distribution is 
given by 



N(Z) 



Njz) 

m 



+ 5N 1 iz) 



(8) 



where / is the Sommer-Larsen factor, and 

SN!(z) =G\iz)-Gliz). (9) 

Gj is the gaussian function adjusted to N(z)/f(z), and G\ 
is its deconvolution, assuming a — 0.20 dcx. Once we find 
N(Z), the SFR history can be recovered by equation (4b). 

3 SIMULATING THE SFR RECOVERY 

We can test the validity of the assumptions of the method by 
means of simulations of the chemical evolution of a galaxy 
experiencing an irregular SFR. We assume a galaxy with no 
infall, and with a prescribed relative SFR, b{t). The initial 
mass function £(logM) and b(t) are connected by the pro- 
cedures described by Miller & Scalo (1979). We have used a 
present-day mass function similar to that of the solar neigh- 
bourhood (Scalo 1986), and the stellar lifetimes of Bahcall 
& Piran (1983). 

The functions £(log M) and b(i) are suitably converted 
into their corresponding quantities according to the formal- 
ism developed by Tinsley (1980), (j>(M) and ip(t), by the 
equations 



M 



0(M) = 
and 



with 



log e ^ (log M) 



M 



M = log e 



JM mi „ 



g M) dM, 



(10) 



(11) 



(12) 



where we take M max = 62M and M min = 0.1M Q (Miller 
& Scalo 1979). 

The equations for the evolution of the gas mass (m g ) 
and metallicity (Z) are given by 



(13) 
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Figure 1. Simulations showing the original (solid lines) and recovered (dashed lines) SFR. 



- [m g (t) Z(t)] = -Z(t)[l - R\m + y[l - R]^{t), (14) 

where y is the metal yield, and R is the mass fraction re- 
turned to the interstellar medium. 

The AMR for this galaxy is directly obtained by solving 
the above equation for Z(t). The metallicity distribution 
is derived by the following procedure. Assuming a sample 
comprising iV t ot stars, the number of stars born at (t, t + At) 



is 

B(t)At = b(t) — At. (15) 
Tg 

The time interval At corresponds to a metallicity interval 
AZ(t). For each star born in (t, t + At), we associate a ran- 
domly chosen metallicity between [Z(t),Z(t) + AZ(t)]. To 
simulate the intrinsic cosmic scatter, this metallicity is fur- 
ther shifted by a random gaussian deviation a = 0.20 dex. 
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We have used several input parametrizations for the 
SFR history corresponding to cases a-h, as shown in Fig. 1. 
Each parametrization is used in the derivation of the age- 
metallicity relation and the metallicity distribution. These 
last two functions are used in order to recover the original 
SFR history. The recovered SFR histories are also shown in 
Fig 1, compared with the original histories. 

Fig. 1 includes several examples, ranging from a con- 
stant SFR (fig. la) to several burst or series of bursts (figs, 
lb-g), and also including epochs of particularly low SFR 
(figs. lg,h). We see from Fig. 1 that there is a good agree- 
ment between the original SFR histories (solid lines) and 
the recovered ones (dashed lines) in the whole time range, 
except at the extremes t = and t = 13 Gyr. Fig. la cor- 
responds to an originally constant SFR. The recovered SFR 
is nearly constant, showing a certain degree of noise. This 
could mask less intense star formation events, making them 
indistinguishable from a constant star formation activity. 

The bursts of Figs. lb,c occur at the same time, but the 
intensity of the burst in Fig. lc is six times lower than in Fig. 
lb, and its duration is almost half the duration of the larger 
burst. While the recovered SFR in Fig. lb does indicate the 
occurrence of the burst, in Fig. lc it is barely noticeable. 
The method is likely to smear out the star formation events, 
in the sense that large amplitude bursts are recovered as 
broad-looking less intense features. As a consequence, we 
expect that the method has a finite resolving power, that 
is, star formation events occurring in a very short timescale 
cannot be recovered. This is illustrated by Figs, le and If. 
The original SFR in Fig. le is composed by two bursts with 
their maxima spaced by 2 Gyr, the second being less intense 
than the first. The recovered SFR combines these bursts 
together in a large star formation event. The second, less 
intense burst is noticeable by a small bump at t ~ 8 Gyr 
(Fig. le, dashed line). Fig. If shows a similar case, where the 
bursts are spaced by 3 Gyr. In this case, we can successfully 
isolate both bursts, suggesting that the resolution of our 
method is about 1.5 Gyr. 

A feature shared by all recovered SFR histories deserves 
more attention. All histories show a burst at t = 13 Gyr. 
This burst is probably an artifice created by the method, im- 
plicit in the counting procedure that transforms the metal- 
licity distribution in a distribution of chemical ages [equation 
(4b)]. The method assumes that all stars with z > Z(Tg) 
were born at Tc, which largely overestimates the real num- 
ber of stars born at To- 

A similar problem occurs in the recovery of b(t) at t — 0. 
All SFR histories, except that in Fig. Id, show a scarcity of 
stars born at t = relative to the original histories. In fact, 
the number of stars born at t — does not include the stars 
with [Fe/H] < —1.2, which are removed from the sample, 
according to the chemical criterion adopted by Rocha-Pinto 
& Maciel (1996), in order to separate halo and disk objects. 
As a conclusion, we expect the present method to produce 
correct results with a resolution of about 1.5 Gyr during the 
whole time span of galactic evolution, except at the extremes 
t = and t = 13 Gyr. 



Figure 2. Age - metallicity relations from the literature. 



4 APPLICATION TO THE LOCAL DISK 
4.1 The age-metallicity relation 

For the recovery of the SFR history in the local disk, we have 
used the metallicity distribution of G dwarfs recently de- 
rived by Rocha-Pinto & Maciel (1996). Unfortunately, there 
is no well-established AMR in the literature. There are many 
determinations of this relation, but they are generally not 
consistent with each other, so that we decided to take into 
account several age-metallicity relations. These include the 
empirical AMR of Twarog (1980), Carlberg et al. (1985), 
Meusinger, Reimann & Stecklum (1991), Edvardsson et al. 
(1993), and the parametrization of the AMR given by Rana 
(1991). All these relations are shown in Fig. 2, where we 
take the present time as 13 Gyr. To allow their use in our 
method, all empirical AMRs were parameterized by an ar- 
bitrary increasing function of time. 

The AMR by Carlberg et al. (1985) is in remarkable 
disagreement with the others concerning the early times of 
galactic evolution, suggesting that the gas was initially very 
rich. There is a generalized suspicion that this relation is 
biased toward old metal-rich stars (Meusinger et al. 1991; 
however, see Strobel 1991 and Marsakov et al. 1990). We 
decided to consider the SFR history derived from this rela- 
tion as less conclusive. 

The AMR by Edvardsson et al. (1993) is a byproduct 
of a comprehensive analysis of spectroscopic abundances of 
189 F and G stars. The curve shown in Fig. 2 corresponds 
to average points in 1 Gyr bins, and generally agrees with 
the AMRs by Twarog (1980, 329 stars) and Meusinger et al. 
(1991, 536 stars), except at times roughly in the range 7-11 
Gyr, when the AMR by Edvardsson et al. (1993) gives lower 
abundances by approximately 0.1 dex. 
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Figure 3. Derived SFR histories in the local disk for each adopted AMR. Also shown are the location of the peaks of bursts A, B and 
C (Majewski 1993), and the region corresponding to the Vaughan-Preston gap. 



From the work of Edvardsson et al. (1993), the cosmic 
scatter in the insterstellar medium varies between 0.10 to 
0.28 dex, with a most probable value of 0.20 dex. This value 
agrees well with a previous estimate by Twarog (1980). In 
our work, we will generally adopt a somewhat lower value, 
g = 0.15 dex. For a given AMR, the method for the recovery 
of the SFR masks almost completely the chemical evolution 
if we assume an abundance scatter of the same order of the 
dispersion in the metallicity distribution, which is ~ 0.20 
dex according to Rocha-Pinto & Maciel (1996). 

4.2 The local SFR history 

In Fig. 3, we show the SFR histories resulting from the use of 
the adopted AMR, again with the present time correspond- 
ing to t — Tq = 13 Gyr. We also show the location of the 
peaks of bursts B and C (Majewski 1993), and the region 
corresponding to the Vaughan-Preston gap. 

Several remarks can be made on the recovered histo- 
ries. First, we will consider the histories derived from the 
AMR relations by Meusinger et al. (1991), Twarog (1980), 
and Rana (1991), which show a better agreement with each 
other. In general, these histories show the occurrence of an 
extended enhanced star formation era 6 to 10 Gyr ago. We 
can also see a "lull" between 1 to 3 Gyr ago, when the star 
formation activity has decreased considerably. The results 
also indicate the occurrence of an intense burst at t = Tq- 
However, from the discussion in the previous section, wc 
expect most of this burst to be an artifice produced by our 
method. It should also be noted that the history from Rana's 
AMR cannot be taken independently from that derived from 
Twarog's AMR, as the parameterization presented by Rana 
(1991) is based on the data by Twarog (1980). 



Figure 4. Comparison of the results using the AMR of Edvards- 
son et al. (1993) and different values for the cosmic abundance 
scatter. 

The large star formation event present in these histories 
at 4 to 5 Gyr (8 to 9 Gyr ago) clearly corresponds to burst C 
(Majewski 1993). There is also an indication of another star 
formation enhancement at 7 to 8 Gyr (5 to 6 Gyr ago) Gyr, 
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corresponding to burst B (Majewski 1993). This feature is 
very similar to the result of our simulations shown in Fig. 
le, where the occurrence of a less intense burst close to the 
main burst is recovered as a small bump. 

The history derived from the AMR by Carlberg et al. 
(1985) shows an essentially flat SFR, apart from strong 
events at t = and t = 13 Gyr. We should recall that the 
AMR of Carlberg et al. (1985) predicts a high initial metal- 
licity for the disk. By analogy with the problem at Tg, the 
method also assumes that all stars with [Fe/H] < [Fe/H] t=() 
were born at t = 0. As long as [Fe/H] t=0 is very high in 
this AMR, it can be understood why this SFR history has 
so many stars at t = 0. In view of the fact that the AMR 
of Carlberg et al. (1985) is probably biased, we decided to 
disregard the corresponding SFR history. 

A more complex situation occurs with the AMR by Ed- 
vardsson et al. (1993). From the results of Fig. 3, the history 
derived from this AMR seems to indicate a nearly constant 
SFR over the disk lifetime, except that some remnants of 
the star formation events shown in the other histories are 
still visible at about 5 and 10 Gyr. A more detailed analysis 
of this SFR history is shown in Fig. 4, where we compare 
the results using two different values for the cosmic scatter, 
namely a = 0.15 and a = 0.20. In this figure, it can be 
seen more clearly that some star formation enhancements 
are present, which approximately correspond to bursts B 
and C (Majewski 1993). Therefore, the results derived from 
the AMR by Edvardsson et al. (1993) do not contradict the 
results derived from the remaining AMRs, namely those by 
Twarog (1980), Meusinger et al. (1991), and Rana (1991), 
although the derived SFR enhancements are considerably 
smaller. 

We obtained an average SFR history using the derived 
histories, excluding that derived from the data of Carlberg 
et al., as shown in Table 1 and Fig. 5. The formal errors are 
shown, and the last temporal bin is given within parentheses, 
owing to the limitations of our method. 

Fig. 5 also shows a comparison of our average SFR his- 
tory and some previous determinations in the literature: 
Twarog (1980) and Barry (1988), as presented by Noh & 
Scalo (1990); Meusinger (1991a), after transforming to rel- 
ative SFR; and Soderblom et al. (1991), after corrections 
by Rana & Basu (1992). Note that the main events found in 
our SFR history are also present in the other determinations. 
Particularly, we note a large star formation era from 3 Gyr 
to 9 Gyr, and a considerably decrease in the star formation 
activity between 10 and 12 Gyr. 

The second burst at 8 Gyr, suggested by Fig. 4, which 
corresponds to Majewski burst B, is not reproduced in the 
SFR mean history, as it is not present in the histories derived 
from Meusinger et al. (1991) and Twarog (1980). However, 
we do not discard the possibility of its occurrence, as long 
as all other histories in Fig. 5 seem support it. It is possible 
that this burst is not so intense as to be resolved by our 
method. 



4.3 Concluding Remarks 

Our results are in very good agreement with some previ- 
ous determinations of the star formation history. However, 
it must be stressed that the mean history we have obtained, 
as well as the other histories presented in Fig. 5, reflect the 



Table 1. Mean SFR history 
for the local disk. 



t (Gyr) 


b{t) 


to 1 


1.14 ± 0.20 


1 to 2 


0.79 ±0.15 


2 to 3 


1.09 ±0.15 


3 to 4 


1.26 ± 0.22 


4 to 5 


1.64 ±0.39 


5 to 6 


1.42 ±0.10 


6 to 7 


1.14 ±0.24 


7 to 8 


0.95 ±0.08 


8 to 9 


0.78 ±0.13 


9 to 10 


0.61 ±0.25 


10 to 11 


0.54 ±0.28 


11 to 12 


0.39 ±0.25 


(12 to 13 


1.28 ±0.22) 



Figure 5. Comparison of our mean SFR history and previous 
determinations in the literature. 

history of the star formation amongst late-type dwarfs. We 
have not made any assumptions regarding possible varia- 
tions of the initial mass function (IMF) during the star for- 
mation bursts. In the case that the IMF favours massive 
stars during the bursts (Mezger 1987), we would expect that 
the global SFR history should be very different from that we 
have recovered. As long as there are few evidences concern- 
ing variations in the IMF in these environments (Scalo 1986, 
1987; Larson 1992), we believe that the history we have re- 
covered is representative of the global SFR activity. 

There are many independent evidences for the occur- 
rence of star formation bursts in the solar neighbourhood. 
As long as the stars presently observed in our vicinity are 
likely to have born at different galactocentric distances, the 
bursts we find among solar neighbourhood stars could have 
occurred over a larger part of the disk (Meusinger 1991b). 

Majewski (1993) gives a scenario for the evolution of 
our Galaxy with these bursts. The author suggests that the 
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burst occurs during a certain time after which the global 
gas heating, due to supernova explosions, inhibits further 
star formation. During an interval of a few Gyr, the hot 
gas is efficiently mixed, so that the next stellar generation 
(characterized by the next burst) has a greater mean metal- 
licity and a smaller abundance scatter. This scenario can 
satisfactorily explain the chemical properties of the stars 
which define bursts A, B and C (Majewski 1993). Majewski 
also presents some similarities between his stellar population 
scheme as discrete burst populations and the classification 
system proposed by Roman (1950, 1952). In this classifica- 
tion, the disk late-type dwarfs are subdivided into a spectral 
sequence composed by the groups weak CN, weak-line and 
strong-line. The existence of these discrete stellar groups is 
more easily explained by the assumption that the building 
of the disk was not a continuous process. In fact, Majewski 
points out that the chemical and kinematical properties of 
the stars in the groups weak CN, weak- line and strong- line 
are in remarkable agreement with the properties of the stars 
which define bursts C, B and A, respectively. 

A previous explanation for the occurrence of bursts in 
the disk was given by Scalo (1987). He remarks that the sim- 
ilarities between the ages of the bursts in the disk and those 
in the Magellanic Clouds may suggest that the bursts could 
be the result of past interactions between the Clouds and our 
Galaxy. In fact, it is well known that close encounters be- 
tween galaxies could trigger star formation bursts (Combes 
1987). This hypothesis also seems very attractive, because 
orbital calculations show that the Clouds are presently near 
the perigalacticon (Lin, Jones & Klemola 1995), which could 
have triggered burst A. 

An usual criticism with respect to the occurrence of star 
formation bursts in the local disk comes from the traditional 
view that our galaxy is a typical spiral galaxy, where the 
SFR has been constant or slightly decreasing over its life- 
time. In fact, the Galaxy is probably not a typical object, 
and should not be considered so. The existence of nearby 
satellite galaxies (the Magellanic Clouds) and the warp in 
the disk at larger galactocentric distances indicate that our 
Galaxy is rather an interacting galaxy, subject to bursts and 
other phenomena that occur in such galaxies. 
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